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ABSTRACT 

Using archival X-ray data, we find that the catalog location of the X-ray binary Scutum X-l (Set X- 
1) is incorrect, and that the correct location is that of the X-ray source AX J183528— 0737, which is 15' 
to the west. Our identification is made on the basis of the 112-s pulse period for this object detected in 
an XMM-Newton observation, as well as spatial coincidence between AX J183528— 0737 and previous 
X-ray observations. Based on the XMM-Newton data and archival RXTE data, we confirm secular 
spin-down over 17 years with period derivative P w 3.9 x 10~ 9 s s _1 , but do not detect a previously 
reported X-ray iron fluorescence line. We identify a bright (K s = 6.55), red (J — K s = 5.51), optical 
and infrared counterpart to AX J183528— 0737 from 2MASS, a number of mid-IR surveys, and deep 
optical observations, which we use to constrain the extinction to and distance of Set X-l. From these 
data, as well as limited near-IR spectroscopy, we conclude that Set X-l is most likely a binary system 
comprised of a late- type giant or supergiant and a neutron star. 

Subject headings: infrared: stars — pulsars — stars: individual (Set X-l) — X-rays: binaries 



1. INTRODUCTION 

By scanning the Galactic plane with an X-ray payload 
on a sounding rocket, Hill et al. (1974) discovered an un- 
usual X-ray source in the constellation Scutum, which 
they named Scutum X-l (hereafter Set X-l). Set X-l 
was subsequently detected using X-ray instruments on 
the Copernicus, Ariel V, HE AO 1, and Ginga satellites 
(Charles, Mason, & Davison 1975; Marshall et al. 1979; 
Reid et al. 1980; Cooke et al. 1984; Koyama et al. 1991; 
Yamauchi & Koyama 1993); in contrast, it was not de- 
tected, at least unambiguously, in surveys done using the 
Uhuru and EXOSAT satellites (Forman et al. 1978; War- 
wick et al. 1988, Fourth Uhuru Catalog). In summary, 
these observations determined that (i) the X-ray flux 
varies over a range as wide as < 0.3 mCrab to 20 mCrab, 
and this variability occurs on a variety of timescales 
(Charles et al. 1975; Marshall et al. 1979; Reid et al. 
1980; Cooke et al. 1984; Warwick et al. 1988; Koyama 
et al. 1991; Yamauchi & Koyama 1993); (ii) the emission 
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has a hard spectrum with a significant low-energy cutoff 
that implies an interstellar absorption column density of 
> 10 23 cm -2 ; (iii) the source is an X-ray pulsar with a 
111-s pulsation period (Makino 1988); and (iv) the X-ray 
spectrum may show line emission at 6.4 keV, presumably 
from iron K fluorescence (Koyama et al. 1991). Based on 
its similarity to other sources at the same Galactic lon- 
gitude, Koyama et al. (1990) concluded that Set X-l is 
likely an X-ray binary in the Scutum arm of the Galaxy 
at about 10 kpc distance. 

The most accurate information on the celestial location 
of Set X-l comes from Reid et al. (1980), who presented 
positions of X-ray sources measured with the modula- 
tion collimator (MC) on the HE AO 1 satellite (HE AO 
1 A3). They derived a grid of diamond-shaped error 
boxes, but concluded that only two of the diamonds 
were consistent with earlier observations (Hill et al. 1974; 
Charles et al. 1975; Marshall et al. 1979), viz., those 
centered at 18 h 34 ro 49?54, -07°37'51'.'3 and 18 h 33 m 46!48, 
-07°38 / 40'/0 (both B1950; the corresponding J2000 po- 
sitions are 18 h 37 m 32!14, -07°35'13'.'7 and 18 h 36 m 29!10, 
— 07°36'06'.'9). Each error diamond was approximately 
±14" in Right Ascension by ±76" in Declination. 

In this paper we present the results of our analysis 
of archival X-ray observations of Set X-l with ASCA, 
XMM-Newton, and RXTE 10 . We find that the source is 
not in either of the two error boxes chosen by Reid et al. 
(1980), but is instead in the next box to the west. The 
presence of a 112-s period pulsation (having evolved from 
the 111-s period that was found previously) and the coin- 
cidence in position with one of the HE AO 1 error boxes 
make us confident of our identification. We describe the 
X-ray analysis in § 2, and identify a near-IR counterpart 
from the Two-Micron All-Sky Survey (2MASS; Skrutskie 

10 After submission of this manuscript, we were made aware of 
a similar analysis performed by another group (see Helfand et al. 
2004 and http://www.astro.columbia.edu/-jules/axjl835.4-0737/). 
Their conclusions are largely the same as our own. 
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Fig. 1.— Our identification of AX J183528-0737 with the X- 
ray binary Set X-l and the bright (K s = 6.5) near-IR counterpart 
2MASS J18352582-0736501. The image is a red-green-blue com- 
posite of 2MASS K B , H, and J. The contours show smoothed 
ASCA GIS emission, the XMM-Newton position is shown by the 
circle with a radius of 3" (expanded from a 1" radius for visibil- 
ity), and we also show one of the HE AO 1 position diamonds, from 
Reid et al. (1980). 

et al. 2006) and mid-IR surveys in § 3. Finally, we give 
our discussion in § 4 and conclusions in § 5. 

2. ARCHIVAL X-RAY DATA ANALYSIS 

We found that the X-ray source AX J183528-0737, 
discovered by Sugizaki et al. (2001) in their survey of the 
Galactic plane with ASCA (this specific observation is 
from 1997 October 13), was close to but not consistent 
with the published positions of Set X-l (Reid et al. 1980). 
However, there were too few counts in the ASCA data 
to perform a reasonably sensitive search for periodicities 
and thereby uniquely identify AX J183528-0737: Sug- 
izaki et al. (2001) measured a 0.7-7 keV count rate of 
26.9 kscc -1 GIS" 1 in an observation with 4 ksec of good 
exposure, so there are only rj 220 counts (summed over 
the two GIS detectors). 

A 17-ks XMM-Newton observation of 
AX J183528-0737 (observation 0203850201) was 
performed on 2004 September 18/MJD 53266.26. There 
is one bright source detected in the EPIC-pn and 
EPIC-MOS images (the source is not detected in the 
RGS data). This source, which is located at (J2000) 
18 h 35 m 25!8, -07°36'50" with statistical position uncer- 
tainties of C/l (and absolute uncertainties of < 1"), is 
fully consistent with the position of AX J183528— 0737 
(see Fig. 1). It has a background-subtracted EPIC-pn 
count rate of 0.390(9) s" 1 in the 0.5-8 keV band, and 
is consistent with a point source given the angular 
resolution of XMM-Newton. Figure 1 also shows that 
the positions of AX J183528— 0737 determined from 
both the ASCA and XMM-Newton observations are fully 
consistent with one of the error boxes of Set X-l that 
was derived from the much earlier HE AO 1 observations 
although it does not happen to be one of the two 
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Fig. 2. — Power-spectra of Set X-l. We show Z\ power spectra 
for EPIC-pn (top), EPIC-MOS1 (middle), and EPIC-MOS2 (bot- 
tom) data in the 2-8 keV band. The best-fit frequency from the 
EPIC-pn data (8.860(6) mHz) is plotted as the vertical dashed line 
in all three panels. 

error boxes selected as prime candidates for the source 
location by Reid et al. (1980). 
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Fig. 3. — Pulse profiles of Set X-l in the 0.5—8 keV band, repeated 
twice for clarity. We show two cycles of the events detected in 
the EPIC-pn (black solid line), EPIC-MOS1 (red dashed line) and 
EPIC-MOS2 (blue dash-dotted line) folded on the best-fit period 
of 112.86 s; representative error bars are shown on the right-most 
points. 



2.1. Timing Analysis 

To see whether or not this source is Set X-l, we 
searched for evidence of the 111-s period detected by 
Makino (1988) and Koyama et al. (1991). We extracted 
events (using XMMSAS release 20050815) from the pn 
and MOS data with PATTERN < 4 (singles and doubles) 
and energies between 2 and 8 keV within a radius of 
500 pixels (25"). We barycentered the event arrival times 
and constructed Z\ power spectra (Buccheri et al. 1983). 
Searching the Z\ power spectrum of the EPIC-pn data 
for periods from 75-210 s (all much longer than the frame 
time of 73 ms), we find a very significant peak (Z\ am- 
plitude of 93.6, which gives a probability of 2 x 10~ 41 
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TABLE 1 
RXTE Observation Summary 



ObsID 


Epoch 


Exposure 


Trans. a 


CR b 


Ampl. c 


Signif. d 




(MJD) 


(sec) 


(%) 




(s- 1 ) 




20143-08-01-00 


50619.835381(14) 


10361 


74 


81.98 


9.5(2) 


222 


20143-08-01-01 


50620.757608(13) 


11579 


74 


82.55 


7.5(2) 


259 


20142-02-01-00 


50688.532073(30) 


7840 


68 


59.44 


1.6(2) 


46 


20142-02-03-00 


50765.147085(46) 


8148 


70 


71.82 


0.7(2) 


19 



a Average collimator transmission of Set X-l. 

b Mean count rate including background. 

c Amplitude of fitted sinusoid. 

d Significance of period detection, in units of a. 

in a single trial) at 8.860 ± 0.006 mHz (112.86 ± 0.08 s; 
see Fig. 2), and the power spectra made from the MOS 
data also have peaks at 112.69 ± 0.16 s (MOS1) and 
112.69 ± 0.19 s (MOS2), where the uncertainties have 
been calculated according to Ransom (2001). There is a 
small peak at a Z\ power of 22 in the EPIC-pn power 
spectrum, but this is likely a sidelobe of the main peak. 
No other strong peaks were found in the EPIC-pn or 
MOS1 spectra, while in the MOS2 spectrum there is a 
peak at a period of 100.6 s which reaches a power about 
1/2 that of the peak at 112.69 s. However, given the 
much lower amplitude of the main MOS2 peak (23.6, 
compared to 93.6 for the pn) the secondary peak is not 
very strong overall (single trial probability of ~ 10 -7 ). 
Based on this, we conclude that a ~ 112 s pulse period 
was detected, and that the 100.6 s period is likely either a 
noise spike or an instrumental effect in the MOS2 data. 
Folding the data on the 112-s period, we see (Fig. 3) 
strong, asymmetric pulsations with a pronounced inter- 
pulse, not unlike the pulse profiles shown in Koyama 
et al. (1991). The pn data have an rms pulsed fraction 
of 29%. 

To examine the variability of the X-ray source dur- 
ing the XMM-Newton exposure and to search for back- 
ground flares, we binned the EPIC event data in 200-s 
bins to construct lightcurves, which we show in Figure 4. 
In all three EPIC detectors the background is high with 
episodes of strong flaring for the first 2/3 of the observa- 
tion. This flaring is likely responsible (through teleme- 
try saturation) for the dips in the source lightcurves at 
w 5000 s, as the unfiltered EPIC-pn count rate exceeds 
the limit of 400 s _1 and therefore the cameras entered 
"counting mode" 11 (note that excluding the saturated 
regions has no effect on the period determination) . How- 
ever, even away from the times of strong flaring in the 
background, the source strength appears to vary by a 
factor of w 2; see the small increases near 9000 s and es- 
pecially the variations in count rate after the background 
dies out near 11000 s. These variations are similar to 
those seen by Koyama et al. (1991). We also find that 
the spectral hardness stays relatively constant while the 
overall count rate varies. 

Set X-l has also been observed on a few occasions by 
the Rossi X-ray Timing Explorer [RXTE) satellite. We 
have analyzed portions of the observations in which the 
Proportional Counter Array (PCA) was steadily pointed 
at the target location; other portions of the observa- 



TABLE 2 
Period measurements for Sct X- 



MJD 


Instrument 


Period 


CUT) 




(sec) 


47479.6 


Ginga/L AC 


111.001(4) 


48142.7 


Ginga/LAC 


111.194(3) 


50619.8 


RXTE/PCA 


112.23(4) 


50620.8 


RXTE/PCA 


112.22(4) 


50688.5 


RXTE/PCA 


112.28(6) 


50765.1 


RXTE/PCA 


112.37(7) 


53266.4 


XMM-Newton/EPIC-pn 


112.86(8) 


53266.4 


XMM-Newton/EPlC-MOSl 


112.69(16) 


53266.4 


XMM-Newton/EPlC-MOS2 


112.69(19) 



Note. — Quantities in parentheses arc 1 a uncer- 
tainties on the last digit. Ginga/LAC data are rom 
Yamauchi & Koyama (1993). 

tions included short pointed observations or scans over 
regions around the target location in vain attempts to 
precisely determine the position of the source. The longer 
pointed observations comprise data from two observa- 
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11 See http: / /xmm.vilspa.csa.cs/cxtcrnal/xmm_uscr_support/documcn|jj^^ 



Fig. 4. — Background-subtracted lightcurves as a function of 
time since the start of the observation (0.5—8 keV band, with 200-s 
bins) made from EPIC-pn, EPIC-MOS1, and EPIC-MOS2 data, as 
labeled. We also show the EPIC-pn background (over all energies) 
lightcurve in the bottom panel; this background lightcurve was 
scaled down by a factor of 25 before subtraction from the EPIC-pn 
lightcurve (the EPIC-MOS background lightcurves were similar). 
The bright flaring near 6000 s was responsible for the decrease in 
the source count-rate at the same times due to telemetry satura- 
'tjscopy (§ 2.2) we only used the low background 
ed interval. 
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TABLE 3 

Spectral fits to EPIC data of Sct X-l 



Parameter Fit Type a 

PL TB PL+Fe PL+Fe 



Nn (10 22 cm- 2 ) 8.1(6) 7.8(5) 8.1 8.1(6) 

T/kT (keV) 1.48(13) 30^ 1.48 1.49(14) 

Norm b 1.6(3) 2.19(8) 1.6 1.6(4) 

EW Pe (eV) c ■ • ■ ■ ■ ■ 20(38) 23(40) 

X 2 /DOF 84.5/118 83.6/118 84.2/120 84.2/117 



F x (10- n crg s" 1 cm~ 2 ) d 1.2 1.2 1.3 1.3 



Note. — Quantities in parentheses on the fit parameters (ex- 
cluding the flux) are 1 a uncertainties in the last digit; quantities 
without uncertainties were held fixed for the fit. 
a Fit types are absorbed power-law with (PL+Fe) and without 
(PL) an iron line at 6.4 keV and thermal bremsstrahlung (TB). In 
the first PL+Fe fit, the PL parameters are held fixed at the best-fit 
values from the PL fit. 

b Either power-law normalization, in units of 
10 — 3 photons s -1 cm -2 keV - 1 at 1 keV, or ther- 
mal bremsstrahlung normalization, in units of 3.02 X 
10 — 18 J dVn e nj /(4-Trd 2 ), where d is the distance to the source, n e 
and nj are the electron and ion number densities, and the integral 
is over the source volume (based on the XSpec bremss model). 
c Equivalent width of a putative unresolved Fe line at 6.4 keV. 
d Model flux in the 0.5-10 keV band, corrected for absorption. 

tions in 1997 June, one in 1997 August, and one in 
1997 November, and have exposure times ranging from 
4 to 12 ks (see Tab. 1). Data from all 5 proportional 
counter units (PCUs) in the energy range 4-10 keV were 
binned into 1 s time bins. The resulting light curves were 
fit with a sinusoid plus 3 harmonics over trial periods be- 
tween 110 and 115 s. In all the observations a clear dip 
in x 2 is seen near 112.3 s, although the amplitude of the 
fitted sinusoid varied over a wide range. For each of the 
four observations, the error on the period was determined 
from the range of trial periods with x 2 < xLn + ^X 2 
where A% 2 was the full range of variation of x found in 
a comparable region of period space with no signal. We 
list the final period measurements in Table 2. 

2.2. Spectral Analysis 

We fit spectra to the EPIC data. After filtering for 
single and double events (as above), we excised the data 
taken when the background was high, i.e., prior to 11 ks 
from the start of the observation (Fig. 4) and were left 
with exposure times of 4.5, 5.3, and 5.3 ksec for the 
EPIC-pn, MOS1, and MOS2 respectively, out of initial 
exposures of 16 ksec for pn and 17 ksec for the MOSs. 

We extracted source events from the same 500-pixel 
radii regions that we used for timing. For the back- 
ground, we selected pn events from a circle 2500 pix- 
els (125") in radius located near but not containing 
AX J183528— 0737, as the source was close to a CCD 
boundary and we could not define a proper background 
annulus that would be on the same CCD. For the MOS 
data, where AX J183528— 0737 was located in the mid- 
dle of a CCD, we used a background annulus extending 
from 500 to 5000 pixels (250") in radius. We binned the 
data so that each bin has > 25 counts, and then we gen- 
erated ancillary response files (ARFs) and redistribution 
matrix files (RMFs) for each instrument and fit the data 
in Sherpa 12 

12 Part of the Chandra Interactive Analysis of Observations 
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Fig. 5.— X-ray spectrum of Sct X-l, from the XMM EPIC 
data. The upper panel shows the best power-law fit to the data 
from EPIC-pn (black points), EPIC-MOS1 (red circles) and EPIC- 
MOS2 (blue squares). The fit parameters are given in Table 3. The 
lower panel shows residuals in units of a. 

The data from all three instruments are jointly well-fit 
by an absorbed power-law; we give basic fit parameters 
in Table 3 and show the fit in Figure 5. Attenuation 
factors were calculated using the phabs model (which 
uses cross sections from Balucinska-Church & McCam- 
mon 1992 and solar abundances). The quality of the 
fit is good, as indicated by a low reduced % 2 (0.7), and 
the residuals do not show any significant systematic de- 
viations. The data could also be fit with a very hot 
(kT w 30 keV) thermal bremsstrahlung model, but the 
implied absorption was quite similar and at these tem- 
peratures the thermal model is essentially a power-law 13 . 

Interestingly, Koyama et al. (1991) found evidence in 
the spectrum of Sct X-l for an iron emission line at 
6.4 keV with an equivalent width of 0.2-0.3 keV. While 
not required to achieve a good fit to our data, we exam- 
ined whether the addition of such a line to the spectral 
model would improve the fit. To test this, we fit the 
EPIC data with two models: the power-law model frozen 
to its best-fit values with the addition of a line at 6.4 keV 
(third column in Tab. 3), and an absorbed power-law plus 
iron line where all parameters are free (fourth column in 
Tab. 3). The iron line was assumed to be unresolved for 
both fits. In neither case did the addition of the iron line 
improve the fit significantly. There was a small reduction 
in x 2 from 84.5 to 84.2, but this is not significant. The 
best-fit values of the power-law parameters in the fourth 
fit did not change appreciably from those in the first fit. 
Overall, we find no evidence for a line at 6.4 keV, and can 
set a 90%-confidence upper limit to the equivalent width 
of 96 eV. This is formally inconsistent with the results 
of Koyama et al. (1991). They discussed the difficulty in 
properly subtracting the strong Galactic ridge emission, 
and this could be the cause of the discrepancy, although 
the source could also be variable as we discuss below. 

3. OPTICAL/INFRARED COUNTERPART 

(CIAO), http:/ /cxc. harvard.edu/ciao/. 

13 More complicated models of thermal X-ray emission from op- 
tically thin hot plasmas, such as those of Raymond & Smith (1977), 
give similar results if the abundances of metals are low. 
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TABLE 4 
Optical and Infrared Fluxes of 
2MASS J18352582-0736501 



Fig. 6. — Near-IR color-magnitude diagram (left) and color-color diagram (right) of the fb 2000 stars within 500" of 
2MASS J18352582-0736501, the proposed infrared counterpart to Set X-l, from 2MASS. 2MASS J18352582-0736501 is the red cir- 
cle, the arrows indicates Ay = 5, and we also plot a giant sequence (luminosity class III) at an extinction Ay = 24 and a distance of 1 kpc, 
taken from Cox (2000); this shows the stellar colors that imply the minimum extinction (§ 4.3). The empty region at the lower right of the 
color-magnitude diagram is from the 2MASS sensitivity limit, as indicated by the dashed line. 

3.1. Archival Data 

We found a potential IR counterpart to 
AX J183528-0737 in the 2MASS Point Source 
Catalog: 2MASS J18352582-0736501, which is < 0'.'2 
away from the XMM-Newton source (see Fig. 1). This 
object is rather bright and red, with K s — 6.55 ± 0.02, 
J — K s — 5.51 ± 0.03 and H - K s — 1.97 ± 0.03. It 
is among the brightest of the sources in the 2MASS 
color-magnitude diagram of stars within 500" (Fig. 6), 
and is considerably redder than most sources of similar 
magnitudes. Just on the basis of the magnitude, the 
association with Set X-l is very probable, as we find 
w 10~ 5 sources in this region with K s < 6.6 per square 
arcsecond, so the chance for a random alignment within 
0'.'2 is a negligible 2 x 10~ 6 . 

We also find 2MASS J18352582-0736501 in other 
near- and mid-infrared catalogs: the Deep Near In- 
frared Survey of the Southern Sky (DENIS; Epchtein 
et al. 1999) database, the Spitzer Space Telescope Galac- 
tic Legacy Infrared Mid-Plane Survey Extraordinaire 
(GLIMPSE; Benjamin et al. 2003) data-base, the Mid- 
course Space Experiment (MSX; Price et al. 2001) 14 
Point Source Catalog (MSX6C; Egan et al. 2003), and 
the IRAS Point Source Catalog (PSC). We summarize 
all of the infrared data on 2MASS J18352582-0736501 
in Table 4. 

In DENIS, the source is DENIS J183525. 8-073650, 
and it is detected with J = 11.82±0.06, K s = 6.33±0.10, 
but is not detected in i-band (where the limiting magni- 
tude is ~ 18.5). These magnitudes are slightly different 
from the 2MASS values, even allowing for color transfor- 
mations (Carpenter 2001), but we note that the i^ s -band 
measurement is near the DENIS saturation limit. 

In the GLIMPSE Archive (which is more complete 
but less reliable than the GLIMPSE Catalog), it is 
listed as SSTGLMA G024.3361+00.0657, with fluxes of 
3.49 ± 0.13, 2.61 ± 0.10, 4.17 ± 0.12, and 3.42 ± 0.20 Jy 
at 3.6, 4.5, 5.8, and 8.0 /Ltm, respectively. The varia- 
tion of the fluxes as a function of wavelength is puzzling, 
since a star would generally have monotonically decreas- 

14 See http://irsa.ipac.caltech.edu/Missions/msx.html. 



Catalog/ 


Band Wavelength 


Flux 


Instrument 




(fim) 


(Jy) 


MagIC 


r' 


0.6 


< 3 x 10~ 7 




i' 


0.8 


1.3(2) x 10~ 6 


DENIS 


J 


1.2 


0.030(2) 




A'.s 


2.2 


2.0(2) 


2MASS 


J 


1.2 


0.024(1) 




H 


1.7 


0.40(1) 




K e 


2.2 


1.60(4) 


GLIMPSE 




3.6 


3.49(13) 






4.5 


2.61(10) 






5.8 


4.17(12) 






8.0 


3.42(20) 


MSX6C 


A 


8.3 


5.1 




C 


12.1 


7.8 




D 


14.7 


6.5 




E 


21.3 


3.9 


IRAS 




12 


10.3 



Note. — Values in parentheses arc 1 a statis- 
tical uncertainties on the last digits, if available, 
but note that the GLIMPSE data, especially from 
the 3.6- and 4.5-/im bands, are likely somewhat 
saturated. See § 3 for additional discussion. 

ing values. This is likely due to non-linearity /partial 
saturation of this source, and indeed we note that 
SSTGLMA G024.3361+00.0657 is not included in the 
GLIMPSE Catalog, the data-quality flag indicates that 
no non-linearity correction was applied (bit 19), and the 
fluxes are above the nominal, albeit conservative 15 , satu- 
ration limits of 0.44, 0.45, 2.9, and 1.6 Jy in all bands; the 
radial profiles in the 3.6- and 4.5-/Ltm images also show 
effects of saturation. Finally, in MSX6C, it is source 
MSX6C G024.3359+00.0656, and in the IRAS PSC it is 
IRAS 18327-0739. 

3.2. Optical Photometry 

We performed photometric observations of Set X-l 
on 2006 August 28 with the Magellan Instant Camera 

15 See the GLIMPSE Data Products Document at 
http: / / data.spitzcr.caltcch.edu/popular /glimpse / 20050415_enhanced_vl /Docun 
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Fig. 7. — MagIC r' (left) and i' (right) images of Set X-l. In each case we plot a 30" box, with North up and East to the left. The 
position of 2MASS J18352582-0736501 is indicated by the tick marks. 



(MagIC) at an f/11 Nasmyth focus of the 6.5 m Baadc 
(Magellan I) telescope at Las Campanas Observatory 
in Chile. MagIC is a 2048x2048 SITc CCD with a 
0'.'069 pixel -1 plate scale and a 142" field of view. Expo- 
sures of 1830 s in the r 1 filter and 3630 s in the i' filter 
were obtained. The conditions were near-photometric, 
with 0"7 seeing in r' and 0'.'6 seeing in i'. We reduced the 
data according to standard procedures in IRAF by sub- 
tracting overscan regions, merging the data from four 
amplifiers, and flatfielding the data with twilight flats. 
We astrometrically calibrated the data by measuring the 
positions of 50 2MASS stars and fitting for the transfor- 
mation using ccmap: the fit was characterized by an rms 
error of Of! 07 in each coordinate. 

For purposes of photometry, we performed five 3-s ob- 
servations of the LI 10-232 standard field (Landolt 1992; 
Stetson 2000) in each filter. We transformed the tab- 
ulated Kron-Cousins R and / magnitudes for 12 of the 
stars in that field using the results of Fukugita et al. 
(1996) to r' and i' magnitudes and determined zero- 
points for the observations. These zeropoints agreed with 
the nominal values for MagIC 16 to within our precision 
(« 0.05 mag). 

As seen in Figure 7, we detect an object at the po- 
sition of 2MASS J18352582-0736501 in the i' image 
but not in the r' image. We estimate that the ob- 
ject has r' > 25.2 (3er limit) and i' = 23.64 ± 0.15, 
and we give the corresponding i'-band flux in Table 4. 
The photometry is consistent with the non-detection of 
2MASS J18352582-0736501 in the DENIS i-band. 

3.3. Near-IR Spectroscopy 

To help determine the spectral type of the counter- 
part to Set X-l, we undertook some limited near-IR 
spectroscopy. The spectra were obtained with the near- 
IR spectrograph NIRSPEC (McLean et al. 1998) on the 
Keck II telescope in low-resolution (R ~ 1400) mode and 
reduced using the standard procedures described by Erb 
et al. (2003). The observations consisted of 4 x 1 s expo- 



sures in the K -band (using the NIRSPEC-6 filter), and 
4 x 2 s exposures in the i?-band (using the NIRSPEC-5 
filter), on the night of 2006 October 3. The resolution 
was 15 A in if-band and 10 A in f/-band. We could not 
obtain standard star data from the same night as the 
observations, but we did rough flux calibrations of the 
.ff-band data using an observation from 2004 Septem- 
ber of the A0V star HD 40335 and of the K-b&nd data 
using an observation from 2006 June of the A2V star 
HD 201941. There were variable clouds during the ob- 
servations, and the calibration data are not of the highest 
quality. Therefore the absolute flux scale is uncertain to 
about a factor of 2, and some smaller-scale deviations 
in the continuum level are also present due to imperfect 
correction for atmospheric transmission. 

We show the spectra in Figure 8. We can immediately 
say that there are no strong emission lines such as one 
might expect from the stellar wind in an X-ray binary. 
Moreover, we do not see strong H I Paschen or Brackett 
series absorption that would indicate an early-type star 
- the type of star that is most commonly associated 
with X-ray binaries of this pulse period (e.g., Bildsten 
et al. 1997; Corbet et al. 1999). Instead, we see strong 
CO absorption and lines from a number of metals. Com- 
paring the lines that we identify with the sequences pre- 
sented in Wallace & Hinkle (1997), Meyer et al. (1998), 
and Rayner et al. (2006, in prep. 17 ), we find the closest 
match is with late-K to early-M stars. Given the poor 
flux calibration we were not able to make a quantita- 
tive classification, but the identification of the strong CO 
bands as well as absorption from neutral metals seems se- 
cure; in particular, the comparable strength of the Na I 
doublet and the Ca I triplet along with the presence of 
Al I in the K-band data seem to indicate that the type 
is not too late. We have not attempted to determine the 
luminosity class of Set X-l, as the detailed measurements 
necessary for that are beyond the tolerance of our cali- 
bration, although the strength of the CO bands argues 
for luminosity class I— III (see Meyer et al. 1998). 



16 See http://occult.mit.edU/instrumcntation/magic/#rpt_txt. 



See http://irtfweb.ifa.hawaii.edu/-spex/spcxlibrary/IRTFlibrary.html. 
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Fig. 8.— Near-IR spectra of 2MASS J18352582-0736501 in the H (left) and K (right) bands. We show the spectrum of Set X-l along 
with the comparison stars HD 35601 (spectral type Ml. 5 lab-lb with an added extinction of Ay = 24; from Rayner et al. 2006, in prep.) 
and the X-ray binary GX 1+4 (X-band only; Chakrabarty et al. 1998) as labeled; the fluxes of both comparison stars have been shifted 
arbitrarily. The main absorption lines in the //-band are the Au = 3 CO bands (Goorvitch 1994), along with Mg I and Al I. In the it-band, 
the main lines are Mg I, Al I, Na I, Fe I, and Ca I, although GX 1+4 also shows a strong, broad H I Br 7 emission line; the very sharp 
drop at 2.29 jim is the beginning of the CO u = 2 — > band. The flux scale is approximate. 



For comparison, in Figure 8 we also plot the spec- 
trum of the Ml. 5 lab-lb star HD 35601 (taken from 
Rayner et al. 2006), reddened with A v = 24 (see § 4.3). 
A visual inspection shows that our choice of comparison 
star is reasonable. In the -ff-band the match is partic- 
ularly good, with the overall slope also agreeing. This 
then tells us that Ay ks 24, and although we have not 
corrected for any intrinsic reddening of the comparison 
star, the intrinsic reddening should only be « 1.5 mag 
based on its membership in the Aur OBI association (see 
Levesque et al. 2005). In the A-band the match is not 
as good, with some curvature present in the NIRSPEC 
data from poor calibration, but the absorption features 
and the overall depth of the CO bands agree reasonably 
well. It is possible, though, that the continuum shape in 
the AT-band is a result of water absorption in a very late 
M (later than M7 or so) star rather than poor calibration 
(see Cushing, Rayner, & Vacca 2005). 

4. DISCUSSION 

The detection of pulsations at approximately the same 
period as that found by Koyama et al. (1991), along with 
the position coincidence between AX J183528— 0737 and 
an alternate HE AO 1 diamond (Fig. 1; it is also more or 
less consistent with the other X-ray positions), secures 
the identification of AX J183528-0737 with Set X-l. 
We now discuss the implications of our measurements 
in more detail. 

4.1. Spin-Period Evolution 

Comparing the spin periods that we measure here to 
the values measured by Yamauchi & Koyama (1993) with 
the Ginga satellite (Table 2), we see evidence for secular 
spin-down over 17 years. Yamauchi & Koyama (1993) 
inferred a spin-down rate of P = 3.3 x 10~ 9 ss _1 , al- 
though this was only based on two measurements and 
could not account for possible torque variations or pe- 
riodic (i.e. orbital) changes. Our period measurement 
fits this trend reasonably well (Fig. 9): a weighted fit 
(which is dominated by the Ginga measurements) gives 
v = -3.1 x 10~ 13 Hz s" 1 (P = 3.9 x 10" 9 s s" 1 ). So 



it appears that the spin-down has largely been steady, 
averaged over 17 years, although there certainly could 
have been some small-scale torque variations visible in 
the departure of the RXTE period measurements from 
the long-term trend. This excludes the possibility that 
the period change in Set X-l was due to orbital mo- 
tion as considered by Yamauchi & Koyama (1993), as 
no companion can produce a Doppler shift of 2 s over 
17 yrs, and instead suggests that the spin-down is due 
to external torques. Overall, this scenario agrees with 
the conclusion of Yamauchi & Koyama (1993) that the 
compact object in Set X-l is a neutron star rather than 
a white dwarf (cf. Lopes de Oliveira et al. 2006), which 
would need a much larger torque than those that have 
been observed to produce the measured P. 

4.2. X-ray Spectrum 

Even if we ignore the presence or absence of the Fe line, 
our best-fit spectral parameters are not the same as those 
estimated by Koyama et al. (1991), who found a slightly 
softer power-law (r = 2.0 vs. 1.5 here) and more absorp- 



tion (A H = (2-4) xlO 2 



8x10 cm 2 here). In 



fact, most of the X-ray observations to date (which have 
not been of the same quality as the XMM-Newton obser- 
vation) have inferred column densities above 10 23 cm~ 2 . 
While this could just be an effect of instrumental cross- 
calibration and inconsistent fitting techniques, the dif- 
ference is quite large. This could be due to variations in 
absorption by matter associated with the Set X-l system 
over the 17 yrs between the observations (see Yamauchi 
& Koyama 1993); in particular, variations in Ah are of- 
ten associated with absorption from a variable stellar 
wind (e.g., White & Swank 1984; Chakrabarty & Roche 
1997), and a decrease in Ah can even be correlated with 
a decrease in the Fe equivalent width. Overall, the flux 
has decreased by about a factor of 4 from the Ginga ob- 
servations to the XMM-Newton observations (we find a 
flux of w 0.4 mCrab), and the variability of the pulsed 
amplitude during the RXTE observations over the course 
of several months (Tab. 1) may be as much as a factor 
of 10. The unabsorbed flux in the XMM-Newton obser- 
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Fig. 9. — Spin-down of Set X-l, from Ginga (Yamauchi & 
Koyama 1993), RXTE and XMM- Newton data from Table 2 (as 
labeled). The line is a weighted fit with z> = —3.1 X 10~ 13 Hz s _1 . 

vations implies a luminosity Lx — 1.4 x 10 33 d£ pc erg s" 1 
in the 0.5-10 keV band, where the distance to Set X-l is 
(ikpc kiloparsecs. 

4.3. Optical/IR Counterpart: Constraints on Extinction 
and Distance 

The very red colors of 2MASS J18352582-0736501 im- 
ply a large extinction, as no stars have intrinsic colors 
nearly that red. The reddest main-sequence star listed 
in Cox (2000) has J - K s ps 1 (M5V), while the reddest 
giant star has J — K s ps 1.2 (M5III). If we assume intrin- 
sic colors of J-K s < 1.2 for 2MASS J18352582-0736501 
and that the emission we see is photospheric (no ex- 
cess from a disk or wind), this then implies Ay > 24 
(hence the giant track shown in Figure 6). Such extinc- 
tion makes 2MASS J18352582-0736501 very faint in the 
optical, which is consistent with our measured i' magni- 
tude and with the upper limits at r' and in other bluer 
bands. This extinction agrees with the slope of the H- 
band continuum, although given our calibration uncer- 
tainties that is not a strong statement. The X-ray ab- 
sorption implied by this extinction, using the relation of 
Predehl & Schmitt (1995), is A H > 4 x 10 22 cm" 2 , which 
is consistent with our spectroscopic result. 

The Galactic extinction model of Drimmcl, Cabrera- 
Lavers, & Lopez-Corredoira (2003) predicts that Ay = 
19 at 7.5 kpc, and Ay = 28 at 10 kpc, although at these 
distances and extinctions, the model is not very well con- 
strained. So based on extinction alone, and assuming 
that the extinction is extrinsic to the sources, we would 
estimate a distance of > 8 kpc for Set X-l, which is con- 
sistent with the assertion of Koyama et al. (1990), that 
Set X-l is likely in a spiral arm at ps 10 kpc. 

However, while we know that 

2MASS J18352582-0736501 is heavily reddened, 
with the limited data that we have it is difficult to 
determine its intrinsic colors and stellar type with any 
precision. A rough fit to the i'JHK s photometry (and 
the r' upper limit) is reasonably consistent with a 
late- type star at Ay < 30. This inference is consistent 
with what we deduce from the ncar-IR spectroscopy. 
In general, for a given stellar type the distance to 



2MASS J18352582-0736501 is 

logio 4pc - 0.2 [(K s , obs - O.llAy) + (V — K) - My]-2, 

where if s ,obs = 6.55, Ax/Ay = 0.11, and (V — K) and 
My are the color and absolute magnitude of the star, 
which we take from Cox (2000). We determine Ay by 

(J - K s ) ohs - (J - K) 
v (0.29-0.11) 

where we observe (J — AT s )obs = 5.51, (J — K)o is the 
intrinsic color of the star, and Aj /Ay = 0.29. If we take 
2MASS J18352582-0736501 to be a late-type supergiant, 
we have (for M0I) (V - K) ps 3.80, (J - K) ps 0.9, and 
My ps —5.6, so we find a distance of ps 4 kpc. This is 
slightly smaller than that implied by the extinction, but 
gives a reasonable X-ray luminosity of ps 2 x 10 34 erg s _1 . 

4.4. Comparison with GX 1+4 

Given its infrared and X-ray characteristics, Set X-l 
resembles the 2-minute X-ray pulsar GX 1+4, which is 
a symbiotic system with an M-giant mass donor (David- 
sen, Malina, & Bowyer 1977; Chakrabarty & Roche 1997; 
Chakrabarty, van Kerkwijk, & Larkin 1998). GX 1+4 
generally is quite similar to Set X-l: it has a simi- 
lar spin-period and showed a prolonged period of rela- 
tively faint (< 0.5-2 mCrab) emission with steady spin- 
down (Chakrabarty et al. 1997), although it had pre- 
vious bright periods associated with spin-up. We note, 
though, that we did not see the strong H I emission lines 
one would expect from a symbiotic star (Fig. 8), but this 
may be a result of variability or viewing geometry (see, 
e.g., Masetti et al. 2006). We also note that GX 1+4 
appears somewhat less luminous in the near-IR: GX 1+4 
has K ps 8.1 with Ay ps 5, which imply Mr ~ —5.6 or 
My ps -0.5 (Hinkle et al. 2006). If the Set X-l system 
contains a star like this it would only be 500 pc distant; 
the observed extinction in the optical/near-IR and low 
energy X-ray absorption would be hard to explain, and 
the X-ray luminosity would be very low. 

5. CONCLUSIONS 

It is possible that Set X-l is a nearby low-mass X-ray 
binary with an evolved giant companion. This would 
then imply a low X-ray luminosity — comparable to 
those implied by Masetti et al. (2006) — and that the X- 
ray and optical/IR absorption are intrinsic to the source, 
perhaps caused by material in the stellar wind (e.g., 
Revnivtsev et al. 2003; Filliatre & Chaty 2004). This 
would be consistent with the possible variations in Ah 
seen between our XMM-Newton data and previous obser- 
vations, but those variations are far from robustly deter- 
mined. Additionally, the reasonable agreement between 
the X-ray and optical/IR column densities (c.f. Filliatre 
& Chaty 2004; Kaplan et al. 2006) suggests that the ab- 
sorbing material is largely distributed along the line of 
sight. 

We therefore believe it very likely that Set X-l is an 
X-ray binary with a giant or supergiant late-type com- 
panion located at a distance of > 4 kpc. In this case, 
the mass donor in Set X-l would be larger than several 
hundred Rq (the companion to GX 1+4 has a radius 
Ps IOO-Rq) and significantly more massive than in the 
previous scenario (i.e., a supergiant of ps 10M q rather 
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than an evolved giant of 1-2M ). With such a large 
companion, even at a distance of 10 kpc the X-ray lu- 
minosity would be very low for steady Roche-lobe over- 
flow accretion, so therefore the accretion must be either 
wind-fed or in an elliptical orbit. With cither compan- 
ion (a giant or supergiant), if we assume that the accre- 
tion is wind-fed, we can set a lower limit in the orbital 
period P rb- For a circular orbit and MOI star (with 
stellar radius _R* = 500i? Q and mass M* = 13M Q ), 
we would need P or b > 2 yr so that it was contained 
within its Roche lobe, while a star like the companion to 
GX 1+4 (R* = lOOi? and M* = 1.2M Q ) would have 
^Porb > 0.8 yr. 

With a large companion, Set X-l is a good candi- 
date for an eclipsing system since the companion will 
occult the neutron star with a significant probability 
(Ray & Chakrabarty 2002), and such eclipsing systems 
offer valuable constraints on the neutron star mass dis- 
tribution (van Kerkwijk, van Paradijs, & Zuiderwijk 
1995). Further X-ray observations to measure the or- 
bit of Set X-l and fully study the spectrum, along with 
optical/IR observations to determine the parameters of 
2MASS J18352582-0736501, should be quite valuable in 
unraveling the nature of the system. 
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